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OWGG - G REAL*8 *#* ; , AL*8 power routine 
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004/ le: OTSPOWGG. EDIT: JCW2004 


LARA RAAAAAARARAAAASALALASAL ALES SEER EEE ELE SE SESE EEE EEE EEE TESS CEES C EEE aa 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE IS cea aero UNDER A LICENSE AND MAY BE ee AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 


COOOooooooooooo 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
AND 0 OT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 

DIGITAL ASSUMES NO RESPONS * 


ITY FOR THE USE OR RELIABILITY OF ITS 


IBILITY 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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“ABSTRACT: 
G double base to G double power. 


Floating overflow can occur 
vctasnbhen ex eneny tat ren can occur if: 
yo ve base 
ase and power is 0 or negative. 


) 
| 
| 
FACILITY: Language support library - user callable 


VERSION: 2 
HISTORY: 


AUTHOR: 
Bob Hanek 8-Mar-83: Version 2 


MODIFIED BY: 
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: -SBTTL HISTORY ; Detailed current edit history 
5 ; Edit history for Version 2 of OTS$POWGG 


- ¢ Change fable INDEX to be local instead of GLOBAL. LEB 24-May-1983 
- nesrugy ten us ing INDEX(Rx) to INDEXCRx] to avoid Linker 


; - Change 

3 arrett. ay 
: 2-004 = Changed the ACMASK to include registers RY and R10. This is needed 
; sot SS entry masks of OTSPOWGG and UVXPOWGG are identical. 


; program section for OTS$ code 


Constants 


61 
83 
0 64 
65 
0 96 
4 ° JCW 17-JUL-84. 
909 $3 F 
00 0 -SBTTL DECLARATIONS 
Boe 71 
000 ie 3 
0000 75 ; INCLUDE FILES: 
44 74 ; 
000 75 
0000 76 ; 
0000 77 ; EXTERNAL SYMBOLS: 
0000 78 ; 
0000 79 
0000 80 -DSABL GBL 
0000 81 -EXTRN MTHSSSIGNAL :; Error signal routine 
0000 Ht -EXTRN MTHSK_UNDEXP ; Undefined exponentiation code 
0000 8 eEXTRN MTHSK_FLOUNDMAT : Floating point underflow code 
Bas Hf: EXTRN MTHSK_FLOOVEMAT ; Floating point overflow code 
0000 86 ; 
0000 Hf 3; MACROS: 
0000 8 ; ; 
444 34 : SSFDEF : Define stack frame symbols 
9000 31 ; EQUATED SYMBOLS: 
94444 0900 38 ; base = 4 ; base input formal - by-value 
0000000c 000 94 exp = 12 : exponent input formal - by-value 
00000008 0000 95 fexp = 8 3 exponent when base is floating 
; | 96 3_ by-value 
000007FC 600 $3 ACMASK = “M<R2, R3, R4, RS, RO, R7, RB, RY, R10> 
° : OWN STORAGE: None 
¢ ; PSECT DECLARATIONS: 
0000 4 ; .PSECT _OTSSCODE PIC,SHR,LONG, EXE ,NOWRT 
1) 
12 
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OT spr oucs G REAL*8 ** G REAL*8 power reuttas 19° SEP=-19 r 38:39 Par /yes Be v04-00 
DECLARATIONS ~SEP-1984 MTHRTL SRC 50TS SPOWGG.MAR; 1 
00000000 00000000 00000000 00004002 : if -OCTA osc AOTC Ae 
5404A 172 LOGTAB: .QUAD a iH 2 5404A 3 1319093633 040E+ 
SBOE slenaee HESE Ue tocrn: aug cagpsgononesegen sg. lsorsessannage 
198 4 A 174 “QUAD = ra 43 AC709405 3 307 oae8s $368 E+0 
0000: " 4h A 175 + QUA D yy 00000000000 3 -00000000000000000E +0 
; 178 LOGLEN = <.=-LOGTAB 
636C74B5 78B0E177 882F7652 71544007 3 18 CO: -OCTA %X636C748578B0E177B82F 765271544007 
C 180 
FFO6BA87 30963064 C 18 EXPTAB: .QUAD “XFF06BA8730963D64 $ -14345715681320842E-12 
1AAQEGEF DBgSSDE$ C 1 § eQUAD “X1AA eet ttaitd: 3 0:397372408131 76138610 
7BD56FB9 B ABTES 8 DO) Os «18 » QUAD ann gor BRecnRaEOs : 0-917 ¢306 017354366E-08 
ODDD703 6B083EDC D8 184 -QUAD “X80DDD7036B8083EDC : 3-19? Fete oA TE TA fet 
CS8FFF82 Oo er ae O2EO0 185 -QUAD *X PEE ecBreD orate A. Bet 223813860 24680E-03 
S9OEFFEFA 363 FB6 O2E8 186 -QUAD “X39EFFEFAZE423FB6 i 0.21660849392498290E-01 
00000000 00000000 O02FO 187 - QUAD 440980000000000000 : 0.00000000000000000E +00 
00000007 8 re 133 EXPLEN = <.-EXPTAB>/ 
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EA © G REAL*8 giving 6-SEP-1984 11:28:13 TURTLE SREDOT SPOueG MAR: 1 
~SBTTL OTSSPOWGG - G REAL*8 to G REAL*8 giving G REAL*8 result 


esgron - 


om 


++ 

; FUNCTIONAL DESCRIPTION: 
OTSS$POWGG - G REAL*8 result = G REAL*8 base ** G REAL*8 exponent 
The G REAL*8 result is given by: 


base exponent result 

= 8 >0 0.0 

= = 0 Undefined Exponentiation 
= <0 Undefined Exponentiation 
<0 any Undefined Exponentiation 
>0 >0 2“(exp * Log2(base)) 

>9 #89 pe Ste 

> 0 <0 2“(exp * Log2(base)) 


CSOoOCCOCCCCOOOOCOCOCOOOOOCOoCOoOOoOO 
POPOPOPOPONOPOPoPoNonononopopnorornopnonry 
TNA 


Mroronon 


Floating Overflow can occur. 
Undefined wig ogg hy ie can occur if: 

1 ase is 0 and exponent is 0 or negative 
2) base is negative 


: CALLING SEQUENCE: 


SS SS QOOTTKOOOCOCOCOCOWO0000n0v0no 


MOPROROROROROROROROROnononononononong | 2 es = 


power.wg.v = OTSSPOWGG (base.rg.v, exponent.rg.v) 


INPUT PARAMETERS: 
; base and exponent parameters are call by value 


; IMPLICIT INPUTS: 
: none 
; OUTPUT PARAMETERS: 
: none 
IMPLICIT OUTPUTS: 
: none 
; FUNCTIONAL VALUE: 
; OTSSPOWGG - G REAL*8 base ** G REAL*8 power 
; SIDE EFFECTS: 
SIGNALS floating overflow 
SIGNALS sven tog underflow if underflow detection is enabled 
SIGNALS a EXP (82 = * UNDEFINED EXPONENTIATION') if 
| 


1) Base is 0 and exponent is 0 or negative 
2) base is negative 


NAA WR OS OOD NAME WIN HO OD NAME WIN 9 ODNAU EWN OOO NAUE WR OOONOUS WN WO 
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REAL*8 power Foyt ioe ys a i ae 
EAL*8 to G REAL* 


giving 6-SEP-19 
ENTRY OTSSPOWGG, ACMASK 


Move x to RO r x = 0 and 
3; EXPONENTIATI 
CVTGH exp(AP), -(SP) 
MOV ase(APS, RO 
BGTR DEF INED 
BLSS UNDEF INED 
TSTw (SP) 


BLEQ UNDEF INED 


a 


wo 

n= 
> 
=< 
~ 
< 
= 

con 


acro V04-00 
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=< 0, return ‘UNDEFINED 


- Jif x <0, 0 
ON’ error condition, otherwise attempt to compute x*ty 


Conv 

RO/R1 <== x 
If x > 0 attempt to compute xtty 
Branch to error code for x < 0 

; Test sign o 


; ot pe mag continues here, this implies that x = 0 and y > 0. Return 


Ktty 


RET 


If oe continues here, this implies that an undefined exponentiation 
- Signal error and return 


; was attempte 
UNDEFINED: 

MOVZWL #*x8000, RO 
CLRL RI 


MOVZBL #MTHSK_UNDEXP, -(SP) 
CALLS #1, G*ATHSSSIGNAL 


RET 
: If presne ing continues here will attempt to co 
3; We begin by det 
DEFINED: 

BICL3 #*XFFFFSOOF, RO, RB 

SUBL #*x4010, R8 

SUBL ° 


following identity: 


We begin by determining a as b*i 
and 32 inclusive. Specifically i 
the magnitude of z is minimized. 


We are now ready Se compute log2(x). 


2 
lLog2(2*k*f) = k + log2(a) + ----=- > 
in(2)/ 
where b = 2°(1/32) and i is between 0 
is chosen by table look-up so that 
Since log2(a) = i/32, we may write 


3; Return 


: RO/R1 <== Reserved operand 
; Put error code on stack 


3; Convert error number to to 32 bit 
condition code and signal error. 


quired since there is no JSB entry. 


: NOTE: Second argument is not re- 
: Return 


=< 


B <== 2*44k 


3: RB <-- 2*4e(b 
> RO/R1 <=- f = 


ias 
2*4e(exponent of x 
2e(fraction field 


This computation is based on the 


1 f-a 
son> 2°(2j+1), where 2 = ----- . 
2j+i f+a 


ert exp to h and put on stack 


; Branch to error condition if y =< 0 


mpute xtty as eps. 
ermining k and f such that x = 2“*k*f, where 1 t<d. 


ed exponent of x) 
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oR Oe OTSSPOUGG = G REALSH tore REALSA giving S-SEP=19RG $1:28:7 FATUNTE eae orseouce mars? 729 
y : lLog2(2*k*f) = k + 1/32 + z2*p(2*2). 
5 ; Note that in order to insure an accurate result, log2(2*k*f) must be computed 
19 ; accurately to 68 bits. This will require some H-format arithmetic. 
18 
15 EVAL_LOG2: 
50 50 S6FD 14 CVTGH RO, RO 3 RO/R3 <== f 
54 50 O07 9 6 15 OTL #7, RO, RS 3; Rotate low fraction bits of f 
54 FFFFFF80 8F A A 16 BICL #OXKFFEFEFRO R4 3 R4 <== index to INDEX table 
54 FCBA CF44 90 0341 317 OVB  INDEXCR4], R4 : RG <o= fj 
58 6448 E 47 18 MOVAW nearer R8 3 RB <m- 2°5e(k + 1/32) 
54 50 FD2F CF44 63FD 034 19 SUBHS A_TABLE R4J, RO, R4 3 R4&/R7 <== f - a (NOTE: result is 
035 20 3 exact, i.e. no roundoff error) 
50 06 O32 1 INCL RO : RO/R3 <== 2ef 
50 54 62FD 8 55 § SUBH R4, RO : RO/R3 <-- f + a 
7E 54 50 67FD 0359 DIVH RO, R4, -(SP) ; SP --> 2 
O38¢ 358 
035E $§ 3: Compute 2°5*2*p(2*2) = z2(cO0 + c2"z*2 + c4%2°4 + 66*2%6) 
035E 73 = 2%(c0 + g(z*2)) 
O35E 28 ; 
035E $2 
O35E 0 
56 6E€ 76FD 035 31 CVTH (SP), R6 3 R6/R7 <== 2 (in G) 
50 56 56 45FD 036 3 MULG R6, R6, RO 3: RO/R1 <== 2*2 
FF22 CF 03 50 SSFD 036 POLYG RO, #LOGLEN-1, LOGTAB 3 RO/R1 <== ge(2*2) 
50 50 S6FD O36E 334 CVTGH RO, RO : RO/R3 <== ge(2*2) 
50 FF39 CF 60FD Hes $32 ADDH2 CO, RO 3 RO/R3 <== cO + q*(z*2) 
50 8E 64FD 0378 36 MULH2 (SP)+, RO : RO/R3 <== 2*5ez%p(2*2) 
037C 38 ; 
037C 39 ; Compute log2(x) = k + 1/16 + z*p(z) 
037C 40 ; 
037C 41 
54 58 6EFD 037C rf CVTLH RB, R4 3 Convert 2°5*(k + 1/32) to H 
50 54 60FD 0380 4 ADDH2 R4, RO 3 RO/R3 <-- 2°5*Llog2(x) 
0384 44 
0384 45 ; 
tt ry 3 We now compute xtty = 2*Cy*log2(x)] by writing ytlog2(x) as 
0384 ‘8 : ytlog2(x) = I + j/32 + g/32, 
0384 50 ; where I is an integer, j is an integer between 0 and 31 inclusive, and 
384 51; g is a fraction in the interval [-172, 1/2). It will be convenient to 
B r 26 3 make an initial check for overflow/underflow at this time 
4 354° 
50 6E 64FD 384 55 MULH (SP), RO 3 RO/R3 <== 2°58 qtogesa? 
54 50 8000 8F AB 88 2$ BICW #718600, RO, R4 : R4 <-- exp field of 2*5*y*log2(x) 
54 4010 8F B1 3 5 CMPW #*x401 R4 : 
48 19 38 LSS EXCEPTION_1 ; 
28 6BFD 9 5 CVTRHL Rp. R8 3 RB <mw= 2°5e(1 + j/32) 
4 oGFD 99 60 CVTLH RB, R4 : 
50 4 62FD HY o SUBH2 R4, RO 3 RO/R3 <== 2°5%g 


—TOMNOOWZ EM KOCH TOMMOOWZ EF KOCH LOMMOOMS ES KCK TOMMOOWOZ Er KCK TO MMI 


ui~ 
Wino 


oo 


Pal 
al 
m 

oo 

~ 


pa Ap te pon Ay ap A ey ag af vw 


an 
NOOO NNOLOWSLS & 


“ 

> 

on 
NOUS 
GOoOesrovo 
cvovoVvVlvlT0m vv 


“ 
oo 


“ 
oo 
ow VMIVMIVMMIMMIOF OW! 


NO __BOOSCONA™ ao 
n 


eu On 


-wr OoO-u 
wo £oOoO 


16 
- G REAL*8 ** G REAL*8 power ro sles 16-SE 
OTSS$POWGG - G REAL* to'G BEALSB 6-SE 


Al 
Al 


COOoOCoOoO 
>> 
pares 


PrP 
—— 


si 
3A1 


Q3AC 


MMMMMNOVOVVIFMCCCOWWW 


VOVTI OOM SO PP PPP POOH SNMOOOCOCOCOCO 709 FS NNT NWwD & 


=1384 01:38:38 aw Macro v04-00 Page 8 


giving 13 CMTHRTL.SRCJOTSPOWGG.MAR; 1 (4) 


We can now compute 
xtty = 2*Cytlog2(x)] = 2°CI + j/32 + g/32] 
= (2*1)*(A*(B+1)] = 2*I1*(A + A*B], where 


; A = 2*°(j/32) is obtained from the A_TABLE and B = 2*(g/32) - 1 is obtained 
; by a Min/Max approximation 


CVTHG RO, RO 3 RO/R1 <-- 2*5%g 
POLYG RO, MEXPLEN-1, EXPTAB ; RO/R1 <== B = 2%(g/32) - 1. (Chebyshev 
; polynomial approx. of degree 5) 

BICL3 #*XFFFFFFEO, R8, R4 3; R4 <-- index into A_TABLE table 

mMOovO A_TABLECR4], R4 3 R4/R7 <== A = 2°(j/32) 

CVTHG RZ, R2 7 R2/R3 <-- A 

MULG2 R2, RO 3; RO <== AtB 

CVTGH RO, RO ; 

ADDH2 R4, RO 3; RO <== AtB + A 

CVTHG RO, RO : RO <--2*C(j + g)/32] 

BICW #*XiF, RB 3 RB <-= 2°58] 

ROTL #-1, RB, RB 3 RB <== 2°48] 

ADDW R8, RO 3 RO <m= 2*CI+(j+g)/32] 

CMPW RO, #*XF ; test for over/under flow 

BLEQ EXCEPTION_2 ; see what exception is if neg or = 0 
RETURN: RET 3 otherwise return result in RO 


: Handlers for software detected over/underflow conditions follow 


EXCEPTION_1: 
TSTF RO ; if big ARG > 0 goto overflow 
BGEQ OVER s handler, otherwise go to 
BRB UNDER 3 underflow handler 
EXCEPTION 2: 
TSTw R8 ; test sign of I; if I <0 
BGEQ OVER ; go to overflow handler 


; Underflow; if user has FU set, signal error. Always return 0.0 
UNDER: CLRQ RO ; RO = result. 
BBC #6, SFSW_SAVE_PSW(FP), 2$ ; 

; has user enabled floating underflow? 
MOVZBL #MTHSK_FLOUNDMAT, -(SP) ; trap code for ftoeeing under f Low. 
; Geavert to MTHS_FLOUNOMAT 
; (32-bit VAX-11 exception code) 
; signal condition 


CALLS #1, G*MTHSSSIGNAL 
R return 


ET 


™~m 
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: Signal floating overflow, return reserved operand, -0.0 


et et et ee ee et et et et ak ek ak ek ed eM ah ah ah ahah ah ah ad ad ab ah ahvahab ab abtah ah ah asvahiahiastah ah ah ahahaha vad ad 
tt a 8 9 9 = 4 HOOD DODOOOO OOOO OOO OO O00 009 09 0000 09 09 09 09 09 SI NIN INI NSNNS SNIOAA AAA OO 
DONA MNE WO ODO NA UE WIN 2 O OD NAME WIN © OO NAMES WN OOO NAMES WWN CO OONOUSW 


: else process for overflow 
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59 #01) «2OF = 79 (0401 
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16-SEP-1984 01:58:59 VAX/VMS Macro v04-00 Page 9 
to G REAL*8 giving 6-SEP-1984 11:28:13 CMTHRTL.SRCJOTSPOWGG.MAR; 1 (4) 
420 OVER: MOVZBL #MTHSK_FLOOVEMAT, -(SP) ; Put overflow code on stack 
és] ASHQ #15, #T, RO ; RO = result = reserved operand -0.0. 
4 ; ; RO will be copied to signal mechanism 
4 ; vector (CHFS$L_MCH_RO/R1) so it can be 
424 ; fixed up by any error handler 
425 CALLS #1, G*MTHSSSIGNAL : opel condition 
: § RET ; return = RO restored from CHFSL_MCH_RO/R1 
428 END 
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symbol table -SEP=-1984 MTHRTL.SRCJOTSPOWGG.MAR; 1 ) 
ACMASK = 7F 
A_TABLE R 02 
BASE = 4 
co R 
DEF INED 0 R 
EVAL_LOG2 R 
EXCEPTION_1 0 EO R 
EXCEPTION 2 6R 0 
P = 000 £ 
EXPLEN = 000 
EXPTAB 000002C0 R Os 
INDEX BOR 44 f 0 
LOGLEN = 
LOGTAB 00000290 R 8 
MTHS$$S1GNAL eererene x 0 
MTHSK_FLOOVEMAT eeeeeeee x 00 
MTHSK-FLOUNDMAT tteeeeee Xx 00 
MTHSK_UNDEXP eeereree x 00 
OTSSPOWGG 000002F8 RG 02 
VER QO0003FD R 02 
RETURN QOOOO3DF R 02 
SF$W_SAVE_PSW = ppoea toe 
UNDEFINED 00000300 R 02 
UNDER QOOOO3EA R 02 
bocce mena ence noees a 
H Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
> ~ 00000000 ( 0.) 00 ¢ O.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
SABSS$ 00000000 (¢ 0.) O1¢ 1.) NOPIC USR CON ABS LCL NOSHR EXE RD WRT NOVEC BYTE 
_OTSSCODE 0000040D ( 1037.) 02 ¢ 2.) PIC USR CON’ REL LCL SHR EXE RD NOWRT NOVEC LONG 
g-ewer ewe ese cee ecesc ease ecec em 
! Performance indicators ! 
hipeneindeeaneiniccidiimaigiannmniaiints 
Phase Page faults CPU Time Elapsed Time 
Initialization 31 00:00:00.10 00:00:01.07 
C processing 107 00:00: 3-28 00:00:04.73 
Pass 127 88: 3: §°§ 0: 8 10-36 
Symbol iable sort 0 0:00:00.04 0: 00-08 
Pass 2 88 00: e: 1.12 0:00:04.7 
Symbol table output 4 00:00:00.04 0:00:00.04 
Psect synopsis output 3 00:00: 3-b) 00:00:00.03 
Cross-reference output B2: 0: +90 BB BD. 54 6 
Assembler run totals 361 0:00:04.31 0:00:21.02 


ive working set Limit was 1050 pages. 

9152 bytes (18 pages) of virtual memory were used to buffer the intermediate code. 

There were 10 pages of symbol table space allocated to hold 53 non-local and 1 local symbols. 
488 source Lines were read in Pass 1, produc tng 14 object records in Pass 2. 

9 pages of virtual memory were used to define § macros. 
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SUT Recre Run Statistics RCO TSPOWGG.MAR; 1 


! Macro Library statistics ! 


deem er mee metro ere ese eco meres} 


Macro Library name Macros defined 
~$255$DUA28: (SYSLIBISTARLET.MLB;2 Seer a 
88 GETS were required to define 4 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/D1SABLE=(GLOBAL , TRACEBACK) /LIS=L1S$:OTSPOWGG/OBJ=O0BJ$:OTSPOWGG MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MSRC 
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